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Abstract 



\ Integration by parts is used to reduce scalar Feynman integrals to master integrals. 

Ph. 

^ ; 1 Introduction 

^ . Solution of any problem in perturbative quantum field theory includes several steps: 
> 

• 1. diagrams generation, classification into topologies, routing momenta; 

^ ■ 2. tensor and Dirac algebra in numerators, reduction to scalar Feynman integrals; 

^ ■ 3. reduction of scalar Feynman integrals to master integrals; 

4. calculation of master integrals. 

For a sufficiently complicated problem, all of them must be completely automated. The 
^ . number of Feynman diagram in a problem can be very large. They can be classified into a 
, moderate number of generic topologies. In each of them, a large number of scalar Feynman 
integrals is required (especially if expansion in some small parameter up to a high degree 
is involved). They are not independent: there are integration-by-parts (IBP) recurrence 
relations. These relations can be used to reduce all these scalar Feynman integrals to a 
small number of master integrals py (see also textbooks [2], |3]). The number of master 
integrals can be proved to be finite [1], but the proof is not constructive — it provides no 
method of reduction. 

IBP relations are sometimes used together with recurrence relations in d [SI E]; we 
don't discuss these relations here. Several methods of calculation of master integrals also 
require the reduction problem to be solved: differential equations ^ i8j (see also Chapter 7 
in [2] and the review [9]), recurrence relations in d [6], gluing [10]. Other methods, e.g.. 
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Mellin-Barnes representation, don't depend on IBP. Here we shall discuss only reduction 
(step 3); methods of calculation of master integrals is a separate (and large) topic. 

If the problem involves small ratios of external momenta and masses, there is an addi- 
tional step — expansion in these ratios using the method of regions [H]. It can produce 
new kinds of denominators in Feynman integrals. It is done after (or before) the step 2. 



2 Feynman graphs and Feynman integrals 

Let's consider a diagram with external momenta pi, . . . , pe- If we are considering a generic 
kinematic configuration, this means that the diagram has E + 1 external leg (it can have a 
larger number of legs, if we are interested in a restricted kinematics with linearly dependent 
momenta of these legs). An L-loop diagram has L loop (integration) momenta ki, . . . , k^. 
Let qi = ki, . . . , fc^, pi, . . . , pe {i E [1, M], M = L + E) he all the momenta. The diagram 
has I internal lines with the momenta li, . . . , If, they are linear combinations of g^. Let 
Sij = <ii ■ Qj {j > i)', Ne = E{E + l)/2 scalar products Sij with i > L are external kinematic 
quantities, and 

iV.^(^ + i£ (2.1) 

scalar products with i < L are integration variables. 

Massless and massive lines in a Feynman graph for a scalar Feynman integral are 



•v/vwwwv/v* 



In HQET, the propagator 



-> -P -iO I m? -P - iO 
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(2.2) 



(2.3) 



appears, where the heavy-particle velocity w (f^ = 1) appears in the Lagrangian, see the 
textbooks [T2l [13]. Similar propagators appear in SCET, but with = 0. In NRQED, 
NRQCD the denominator is more complicated: —2MI ■ v + {I ■ f )^ — P — iO. Instead of 
using effective field theories, one may follow a more diagrammatic approach and expand 
the ordinary propagators in specific regions of ki [11]. In all cases, the denominator of a 
propagator is quadratic or linear in the line momentum /. 

The choice of the integration momenta ki is not unique. Momentum conservation at 
each vertex is the only restriction. Some examples of different momentum routings for 
a single diagram are shown in Fig. [TJ It is not easy for a program to recognize that two 
Feynman integrals can be transformed into each other by linear substitutions of ki. In fact, 
the freedom of choice is much larger than suggested by Fig. [1], because all the momenta p, 
ki, ^2 can flow along all lines with some continuous weights. The value of a given Feynman 
integral cannot depend on a choice of momentum routing (all choices agree with Feynman 
rules of the theory, and are equally good). 

Many diagrams have symmetries. For example, the non-planar self-energy diagram 
in Fig. [2^ may be reflected in a horizontal mirror. It also may be reflected in a vertical 
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ki ki — p 




Figure 1: Momentum routings 

mirror; this operation changes the external momentum p — )■ —p, and hence, in order to get 
an identical integrand, we have to do the substitution ki ^ —ki of the integration momenta. 
One more symmetry is less obvious: you can keep the 3 left vertices at rest, and rotate the 
3 right ones around the horizontal axis by vr. The diagram in Fig. |2b (where all the lines 
have the same mass) is very symmetrical. It can be represented as a tetrahedron, and thus 
has the full tetrahedron symmetry group: rotations by mr/S around 4 axes and reflections 
in 6 planes. It is difficult for a program to find symmetries of a Feynman integral, because 
symmetry transformations often have to be followed by integration-momenta substitutions 
to write the transformed integral in its original form. 




a b 



Figure 2: Symmetric diagrams 

We shall consider Feynman graphs containing only 3-legged vertices (generic topolo- 
gies). They have the maximum number of internal lines I. All the remaining diagrams 
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can be obtained from these ones by shrinking some internal hnes, i.e., raising some de- 
nominators to the power 0. They belong to reduced topologies. If the diagram we want 
contains a 4-legged vertex, we can split this vertex into 2 3-legged ones with an internal line 
between them, and say that this internal propagator is shrunk; thus, our diagram becomes 
a particular case of a generic one with power of some internal line. This generalization 
is not unique: a diagram belonging to a reduced topology can be obtained from several 
different generic topologies by shrinking some internal lines. 

For {E + l)-legged tree diagrams, the number of internal lines is I = E — 2 (just 1 
vertex has E = 2 and / = 0; adding an external leg increases / by 1). Adding a loop 
(i.e. connecting 2 points on some hnes) increases / by 3. So, / for {E + l)-legged L-loop 
diagrams is 

Therefore, if L > max(2, 5 — 2E), then there are more scalar products than denominators 
of propagators. In such a case, it is not possible to express all Sij as linear combinations 
of the denominators. For self-energy diagrams {E — 1), this happens starting from L — 3 
loops; for diagrams with more legs — starting from L = 2 loops. 

Vacuum diagrams {E = 0) have to be considered separately. The simplest diagram 
with 2 3-legged vertices has L = 2 and / = 3. Adding a loop increases / by 3, and hence 

/ = 3(L-1), iv_/= (2.5) 

Scalar products which cannot be expressed via denominators appear starting from L = 4 
loops. 

We want all scalar products Sij with i < L to be expressible as linear functions of the 
denominators Da- Therefore, we add irreducible numerators (linear functions of Sij) -D7+1, 

L M 

Da = J2T.^aS^J+ml (2.6) 

i=l j=i 

The only requirement is that the full set Di, . . . , Dj^ is linearly independent. Then we 
can solve for Sif 

N 

^ij-T.'^UDa-ml). (2.7) 

a=l 

If we view the pair (ij) with i G [l,L] and j > i as a single index (it has different 
values), then the matrix A^j is the inverse of A^J. 
Now we define the scalar Feynman integral 

I{ni,...,nN) ^ jr^j^ I d'^ki---d'^kLf{ki,...,kL,Pi,...,PE), 

^'^^ > ^ (2.8) 

f{ki, ...,kL,Pl,...,PE) = -piHi fy^ ■ 



D1^--- D 



N 
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Of course, for irreducible numerators Ua < {a & [I + 1, N]) . The argument of / is a point 
in A'^-dimensional integer space. 

If a diagram contains a self-energy insertion into some internal line, then 2 lines carry 
the same momentum /, and can be joined: 




These 3 different graphs correspond to 1 Feynman integral. 

If the lines around a self-energy insertion have different mass (e. g., 7 and Z^), e.g., 




then their denominators are different but linear-dependent. We shall see that it's easy to 
kill one of these lines using partial-fraction decomposition. It's best to do this before further 
calculations; otherwise, expressing scalar products in the numerator via the denominators 
becomes non-unique. However, this is not possible if these hues are raised to non-integer 
powers (see below). 

If some subdiagram is connected to the rest of the diagram only at 1 vertex (and has 
no external legs) , it can be completely separated and considered a vacuum diagram: 




(2.10) 



In particular. 




(2.11) 



If such a subdiagram has external legs, it still can be separated. The momentum enters 
the rest of the diagram through a new external leg: 




(2.12) 
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In particular, 




(2.13) 



Suppose some subdiagram is connected to the rest of the diagram only at 2 vertices, 
and has no external legs (let's call this subdiagram a). We can always choose the total 
momentum flowing through this subdiagram to be ki. Let the loop momenta k2, . . . , live 
in the subdiagram a; then the remaining loop momenta /ca+I; • • • , live in the remaining 
subdiagram b. Suppose the total degree of the numerator in mixed scalar products ki ■ qj 
(with i < A, j > A) is n. Then the subdiagram a is an integral in /c2, . . . , with n tensor 
indices, depending only on ki. If we detach this subdiagram and then attach it the other 
way round, this means ki — )■ —ki. Let's make the substitution ^2 — ^ —k2, • • • , — ^ —kA- 
All the denominators in the subdiagram a remain unchanged, and 




(2.14) 



The second equality in (12. 9 p is a particular case of this property. 

The A^-dimensional integer space in which Feynman integrals (12. 8 p live can be subdi- 
vided into sectors (Fig. [3]). If Ua > 0, Da is in the denominator; if Ua < 0, it is in the 
numerator (the line is shrunk). The set of sectors is partially ordered: the ++ sector is 

higher than H — which is higher than ; the ++ sector is higher than — h which is higher 

than ; but the sectors H — and — h cannot be compared (Fig. [3]). Each sector has a 

corner — the point with = 1 (if ra^ > in the sector) or ria = (if < in the sector). 
Generally speaking, there are 2^ sectors. But for irreducible numerator D^, we always 
have Ua < 0, and opposite sectors don't exist. Some sectors are trivial, i.e. / = in all 
points. At least, the pure negative sector (where all ria < 0) is trivial. Often there are more 
trivial sectors, when shrinking some lines produces a scale-free vacuum subdiagram (12. lip 
(some sectors are trivial only at some speciflc kinematics, e.g., ( I2.13P ). In sectors just 
above trivial ones (i. e. when the diagram vanishes after contracting any line) a general 
expression for the integral (usually via F functions) can often be obtained. 

Some sectors are transformed into each other by symmetries. All the denominators of 
the original integral become those of the symmetric integral, possibly, after an appropriate 
integration-momenta substitution. The numerators of the original integral (i. e.. Da whose 
ria < 0), after the symmetry transformtion and the loop-momenta substitution, can be 
expressed as linear combinations of the numerators Da in the new sector. Expanding the 
product of powers of such numerators, we can express the original integral as a linear 
combination of integrals in the new sector having different powers of the new numerators. 
A sector may be transformed into itself by some symmetries. 
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rii 



Figure 3: Sectors in the A^- dimensional integer space 

Let's consider a self-energy insertion into a massless line ( 12 ■2p (the HQET propaga- 
tor (12.31) is also "massless" in this sense, because it contains no dimensional parameters 
except /). If this insertion contains no massive lines, then, by dimensionality, this insertion 
just shifts the power of the propagator: 




(2.15) 



where L is the number of loops in the insertion, and the integer n is determined by the 
powers of the denominators in it. A diagram with such insertion [s] can be considered as a 
lower-loop diagram (12.81) with non-integer (rf- dependent) power [s] of some denominator [s]. 
Such a power cannot be compared to integers, and hence there is no subdivision into sectors 
along the corresponding n axis. This factor is always in the denominator, and can newer 
become a numerator. I know no circumstances when a massive propagator gets raised to 
a non-integer powei0. 

^Except the situation when the d-dimensional space is separated into an n-dimensional subspace and a 
transverse [d — n)-diniensional one; integrating a massive propagator in transverse momentum components 
produces a massive propagator in the rt-dimensional subspace raised to a non-integer power. 
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3 Integration by parts 

The Feynman integral f l2.8p does not change if we do a substitution 

ki y MijQj = Aijkj + BijPj 

of its integration momenta, where 

/An ■■■ Bu ■■■ B,E \ 

M= : •.. : : •.. : 

y Ali ■ ■ ■ All Bli ■ ■ ■ Ble j 
is an L X M matrix, provided that the substitution is invertible: 

det A ^ . 

Such substitutions form a Lie symmetry group of the Feynman integral 

Let's consider an infinitesimal transformation fcj — )■ fcj + aqj. The integrand 

as 

f ^ f + aqj- dif . 
If j = 2, also the integration measure changes: 

d% (1 + ad)d% . 

These infinitesimal transformations form the Lie algebra [1 



d'^ki-'-d'^kL OiJ = 0, 



d 



Oij = di ■ qj {i < L, j > i) , di = — . 

dqt 

The generators obey the commutation relation 

[Oij^Oi'ji] = 6ij'0i'j — 6i'jOij' 

(structure constants). 

Let's find an explicit form of the operator 



M 



Oij = d6ij + qj ■ di = d6ij + + Smi)smj 



d 



m=l 



dSr. 



N 



9 ^ A"^' ^ 



a=l 
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(the extra factor 2 at m = i comes from diSu = 2ki). We assume that whenever Sij with 
i > j appears in an equation, it is immediately replaced by Sji] the same holds for d/dsij, 
Afj, Ai^. lij<L {qj = kj), then from we have 

N N M r. 

O^^ = d5,, + + - ^')^ 5 (3.8) 

a=l b=l m=l " 

if j > L {qj = Pj-i), then Smj with m,j > L are external kinematic quantities: 



N 



a=l 



,m=l b=l m=L+l 



' (3.9) 
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The operator 8/ dDa acting on the integrand / fl2.8p raises the power Ua by 1 (and multiplied 
the integrand by lowers n;, by 1. 

Let's introduce operators which act on functions of integer variables and produce 
new functions: 

(n„F)(ni, ...,?!„,..., Um) = UaF^rii, . . . , . . . , tiat) , 

(a+F)(ni, ...,na,.. .^Un) = F{ni, . . . , + 1, . . . , nAr) , (3.10) 

(a~F)(rai, ...,na,.. .,71^) = F{ni, . . . ,72^ - 1, . . .,71^) ■ 

The shift operators are inverse to each other: 

a+a" = a"a+ = l. (3.11) 

They don't commute with the number operators Ua'. 

[a^,n,] = ±5ab^^. (3.12) 

Some authors prefer to use the operators a+ = naa+, 

(a+F)(ni, . . . , Ua, . . . , Un) = naF{ni, . . . , + 1, . . . , Uat) (3.13) 



instead of a+, because 



d 

-a+ , Db^h- (3.14) 



are the only combinations appearing in (13. Sp . (13. 9p . These operators obey the commutation 
relation 

[k-^ ,h-] = Sab ■ (3.15) 

Then the operators Ua = a+a" are not independent; it is sufficient to use a+ and a~. 
Now the identities (13. 5p can be rewritten as the IBP recurrence relations 

0,,(a+,a-)/(ni,...,n^) = -^-^ J d% ■ ■ ■ d'^kLO^jf = , (3.16) 



where the operators Oij are obtained from fl3.8p . fl3.9p using the substitutions (13.141) : 

N N M 

O,, = dS,, - E E E ^r< (l + '5™)a-^ (b- - ml) < L) , (3.17) 



O, 



N 

E 

a=l 



a=l b=l m=l 
N L 



M 



6=1 m=l 



m=L+l 



(z>L) 
(3.18) 



(note the order of and b^). The operators Ojj(a^,a~) obey the commutation rela- 
tions (13.61) by virtue of (13.151) . 

The simplest example is the 1-loop vacuum diagram (Fig. Hj) 



J B"" ^ ' 



D = m^ -k^ -iO. 



(3.19) 



We may set m = 1; the power of m can be reconstructed by dimensionality. The sector 
< is trivial (Fig. Hj). The IBP relation 



{d-2n + 2nl+)V{n) = 



(3.20) 



relate 2 neighbouring values of n (except the n = relation where ^(1) does not appear, 
and we obtain the known fact V^(0) = 0). In the positive sector, we can express V{n) via 
V{n — 1), then V{n — 2), and so on, until we reach V{1). The explicit solution of (13.201) is 



V{n) 



1 r(n-f) 

Tin) r (1 - f ) 



V{1) 



(3.21) 




n 



Figure 4: The 1-loop vacuum diagram and its sectors 

The 1-loop vacuum diagram with masses m and (Fig. [5]) contains linearly dependent 
denominators 

L)i = 1 - fc^ D2 = -k^ , Di- D2 = l (3.22) 
(we set m = 1). In the positive sector, we solve the recurrence relation 

(1-1^ + 2^)7 = (3.23) 
10 



for J; each reduction step decreases rii + n2 by 1, and eventually we get ni = (where 
/ = 0) or ^2 = 0. In the sector n2 < 0, we solve for 2~J; each step decreases \n2\ by 
1, and we end up at r;,2 = or ni = (Fig. |5]). This reduction is equivalent to partial 
fraction decomposition. So, this problem reduces to the previous one (Fig. Hj) (unless n2 is 
not integer). 

^2A 




Figure 5: The 1-loop vacuum diagram with masses m, and its sectors 

Linear dependent denominators often appear in HQET diagrams [15j. For example, 
the vertex diagram in Fig. |6] has 

Di = -2{k + p{)-v = -2{k -v + LOi), D2 = -2{k + P2) ■ v = -2{k ■ v + UJ2) , 

Di - D2 + 2{uJi - L02) = 0. 

The recurrence relation 

[2{ui - U2) + 1" - 2-] / = (3.24) 

allows one to kill one of the heavy lines, similarly to Fig. |5l If an L-loop diagram contains 
more than L HQET lines (with the same v), then its HQET denominators are linearly 
dependent: there are only L scalar products ki ■ v (see, e. g.. Fig. [H]). Some of these HQET 
lines can be easily killed by partial fraction decomposition, unless their powers are non- 
integer. 

The 1-loop massless self-energy diagram (Fig. [7]) has 

D, = -{k+p)\ D2 = -k\ 

It has only one non-trivial sector, and is symmetric with respect to 1 2. We may put 
= —1 (its power can be restored by dimensionality). The d ■ k IBP relation is 

[rf-ni-2n2 + nil+(l-2")] G = 0. (3.25) 
11 
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Figure 6: HQET diagrams with linearly dependent denominators 



If ni > 1, we can use f l3.25p to lower rii + n2 by one; if n2 > 1, the mirror-symmetric 
d ■ {k — p) relation can be used instead. All integrals reduce to 1 master integral G(l, ij^. 




^2 



o o o o 



• • • 



Til 



Figure 7: The 1-loop massless self-energy diagram and its sectors 

In general, when solving the reduction problem in some sector, integrals from lower 
sectors (where some lines are contracted) are considered trivial, i.e., reduction problems in 
those sectors are assumed to be solved. In other words, the reduction problem is recursive: 
to solve it for the fully positive sector, one has to solve it for all the sectors just below it, 
and so on, until we reach trivial sectors (where / = 0). 

Automatic identification of trivial sectors is a useful feature for programs of IBP reduc- 
tion. If there exists a substitution M{a) (13. ip such than the integral is multiplied by a" 
(where 7^ is (i-dependent), then the integral obviously vanishes (contains a scale-free 
subdiagram). The same holds if any numerator is inserted into this integral. Therefore, if 
the integral at the corner of the sector vanishes, the whole sector is trivial. This gives the 
criterion |T3]: if solving all IBP relations at the corner point results in J = 0, the sector 

^Reduction for the 1-loop massless triangle diagram has been considered in [IB] . 



12 



is trivial. This criterion finds all cases like (12. lip , but not sectors which are trivial only in 
specific kinematics, like fl2.13p . 

When considering some specific sector, we should express more complicated integrals 
via simpler ones. To do this systematically, we need to accept some convention which 
integrals are more complicated and which are more simple. This means introducing a total 
order of integer points in the sector. It is natural to require this order to be translation 
invariant: if ni < n2 then Ui + m < n2 + rh. The choice of such an admissible order is, of 
course, not unique. We want the reduction process to move integrals closer to the corner 
of the sector. Therefore, it is reasonable to compare the sums n = Yla l^^l fi^'st: if the first 
integral has a higher n than the second one, the first integral is more complicated. After 
that, we should decide how to order integrals with equal values of n. This can be done in 
different ways (e.g., lexicographically). 

Suppose we have fixed some admissible order in the sector. Let's peek some IBP relation 
and solve it for the most complicated integral. This solution can be used for reduction of 
all integrals in the sector. They either move into lower sectors (when some = 1 and a~ 
is applied), or end up in a hyperplane of dimensionality — 1 in our sector. Indeed, the 
most complicated integral in our IBP relation is naa+J for some a; we cannot reduce Ua 
from 1 to using this relation. So, we are left with the hyperplane = 1. 

Next we update the remaining IBP relations to contain only integrals on this hyperplane 
(and trivial integrals from lower sectors: a~ acting on an integral with = 1 makes it 
trivial). Whenever an IBP relation contains = 2, we lower it back to = 1 using the 
first selected relation. Now the problem is much simpler: the dimensionality of space is 

— 1 instead of A^. Let's select one of those (A^ — l)-dimensional recurrence relations, 
and solve it for the most complicated integral. Then we can use this solution to further 
reduce all integrals on the hyperplane. Here is a catch, however: the coefficient of this most 
complicated integral may vanish on some complicated subset of points, and the reduction 
fails on this subset (this can never happen during the first stage: there this subset is always 
the hyperplane Ua = 1). This subset is still much smaller than the whole hyperplane: it 
consists of some lower- dimensional part[s]. It is better to choose a relation for which this 
subset is simple (e.g., a coordinate hyperplane). 

Next we treat these (A^ — 2)-dimensional subsets in a similar way, and so on. This 
is a sketch of what people usually do when constructing a reduction algorithm by hand. 
This strategy is implemented in a Mathematica program by R.N. Lee based on [Ti] . 
This program is not guaranteed to construct a reduction algorithm for a given topology. 
However, it works successfully (and efficiently) for a large number of highly non-trivial 
examples. It is not publicly available. 

In some rare cases, approaches based on sectors can fail to detect a relation between 2 
integrals in some sector, and declare both to be masters, when in fact they are dependent, 
but this dependency can only be obtained via a higher sector. 
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4 Homogeneity and Lorentz- invar iance relations 



Scalar Feynman integrals obey some relations which can be derived independently of IBP 
but are not really independent — they appear to be linear combinations of the IBP rela- 
tions. 

By dimensionality, an L-loop integral / is a homogeneous function of external momenta 
Pi and masses of degree Ld—2 ^ rii (we assume that all the denominators are quadratic). 
Therefore, 

On the other hand, the derivatives in the left-hand side can be calculated explicitly. Equa- 
tion these two expressions, we obtain a homogeneity relation [1]. 

Let's move all the terms in (14. ip to the left-hand side, and apply this operator to the 
integrand / instead of the integral /. Adding and subtracting the sum over the loop 
momenta fcj we get 

\ i i i / 

\ i i i i / 

The integrand / is a homogeneous function of the momenta qi (both loop and external) 
and the masses of degree —2 ^nj, and this expression simplifies to 




Thus the homogeneity relation is a linear combination of the IBP relations Yl^i ■ ki. 

A scalar integral I does not change if we rotate the external momenta pi. In other 
words, a Lorentz-transformation generator applied to / gives 0. If there are at least 2 
external momenta {E > 2), we can contract this tensor equation with pfp'^ {i ^ j) to 
obtain a scalar relation 




(square brackets mean antisymmetrization) . On the other hand, the derivatives can be 
calculated explicitly. This gives Lorentz- invariance relations p7f]. 

Let's apply this operator to the integrand /. Adding and subtracting the sum over the 
loop momenta ki we get 
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The integrand / is a scalar function of the momenta g^, and hence the first operator gives 
when acting on /. Writing the antisymmetrization exphcitly and commuting derivatives 
to the left, we obtain 



-Pi-knPrT^] f 



9 (4-4) 

= ^-Q^- (PiPj ■ - PjPi ■ kn) f 
n ^ 

(the extra terms from commutation cancel). Thus Lorentz-invariance relations are linear 
combinations of IBP relations |l4j . 



5 Massless self-energy diagrams 

Let's consider the 2- loop integral [ijl (Fig. [8]) 

f{ki, k2,p) = jjn,\jjn, , (5.1) 

D, = -{k,+p)\ D2 = -ik2+p)\ D, = -kl, D^ = -kl, 

D^ = -{h-k2f. 

It is symmetric with respect to the interchanges (1 ^ 2,3 ^ 4) and (1 3, 2 4). 
It vanishes when indices of two adjacent lines are non-positive integers, because then it 
contains a no-scale subdiagram. 




Figure 8: Two-loop massless self-energy diagram (all lines are massless) 

When one of the indices is zero, the problem becomes trivial. If ns = 0, it is the product 

■^The IBP and homogeneity relations for this integral had been derived in [18] slightly earlier than in [1]; 
however, the reduction algorithm had not been formulated. 
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of two one-loop diagrams: 



^3 




(5.2) 



ni 



n2 



If ni = 0, we first calculate the inner loop (it shifts the power of the upper propagator in 
the outer loop), and then the outer one: 



^4 + ^3 + ns - d/2 







f \ 




) 













n2 



(5.3) 



^2 



The cases 77,2 = 0, ns = 0, n4 = are symmetric. If all rii are integer, then all integrals ( 15.2 
are proportional to and all integrals (15. 3p to G2- 




G\ , 




Go 



(5.4) 



where 



Gr 




r (n + 1 -nf) r"+i (f - 1) 

r((n + i) (f-i)) 



(5.5) 



is the n-loop massless sunset integral. 

Now we shall consider the positive sector (all rij > 0). The 82 ■ k2 IBP relation is 

[d-n2-n5- 2^4 + n22+(l - 4") + n^b^iS' - 4-)]G = ; (5.6) 

the di ■ ki relation is mirror-symmetric. The 82 ■ (^2 — ki) IBP relation is 

[d-n2-ni- 2n5 + n22+(l" - 5") + n44+(3" - 5")]^ = . (5.7) 

Let's express G{ni, . . . ,77,5) (with unshifted indices) from this last expression. Each appli- 
cation of this relation reduces rii -|- 77,3 -|- 77,5 by 1 (Fig. [9]). Therefore, sooner or later one of 
the indices ni, ns, 775 will vanish, and we'll get a trivial case ( 15. 2p . (15. 3p . or symmetric to 
it. This means that all integrals in this sector reduce to 2 master integrals (15.41) . Analyses 
of the remaining non-zero sectors is simple but somewhat lengthy; all integrals in these 
sectors reduce to (15.41) too. 
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Figure 9: One step of IBP reduction: projection onto the ni, ns, subspace 



Applying p ■ (d/dp) to fl5.ll) we get the homogeneity relation 

[2{d - - - ^5) -ni-n2 + r2il+(l - 3") + n22+(l - 4")]G' = . (5.8) 

It is nothing but the sum of the 82 ■ k2 relation (15.61) and its mirror-symmetric di ■ ki 
relation. 

Another interesting relation is obtained by inserting (ki+p)^ into the integrand of (15. ip 
and taking derivative d/dp^ of the integral. On the one hand, the vector integral must be 
proportional to p^, and we can make the substitution 

in the integrand. Taking d/dp^ of this vector integral produces fl5.ip with 

inserted into the integrand. On the other hand, explicit differentiation in p gives 

d + ^2(A;i + pf + ^2{k2 + p) ■ {k, + p) , 
jJi U2 

2{k2+p)-{ki+p) = D,-Di-D2. 
Therefore, we obtain the Larin's relation 



|ci + - 723 - 724 - 725 + - ^ 77^) (1" - 3") + 7722+(l- - 5") G = (5.9) 

(three more relations follow from the symmetries). This relation can be used instead of (15. 7p 
to reduce all integrals in the positive sector to the master integrals (15. 4p . see Fig. [9l It is 
surely some linear combination of the IBP relations, but I have no explicit proof of this 
fact. 
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The IBP relations (15. 6p . (15. 7p are particular cases of the triangle relation [T]. Suppose 
a diagram we are considering contains a subdiagram in Fig. (TO) 

Di = ml- {ki + , D2 = ml- kf , 

D3 = ml- {k2 + k3f , Di = ml-kl, D5 = -kl 

(any number of lines of any kinds may be attached to the left vertex; however, it is essential 
that only one line is attached to each of the two right vertices of the triangle). The ^3 ■ k^ 
IBP relation is 

[d-ni-ns- 2n5 + nil+(2- - 5") + n33+(4- - 5")] 1 = (5.10) 

(e. g., (15.71) ). A single application of the triangle relation (I5.10p reduces n2 + + hj 1, 
and this always allows one to kill one of the lines 2, 4, 5. If the line 2 is external instead 
of internal, the lowering operator 2~ becomes just the factor {ml — k\) (where ki is now 
an external momentum), and there is no subset of n's whose sum is reduced (e. g., (15. 6p ). 




Figure 10: Triangle subdiagram 



There are 3 generic topologies of 3-loop massless propagator diagrams: 




(5.11) 

Each has 8 denominators. There are 9 scalar products of 3 loop momenta ki and the 
external momentum p. Therefore, for each topology, all scalar products in the numerator 
can be expressed via the denominators and one selected scalar product. IBP recurrence 
relations for these diagrams have been investigated in [Tj. They can be used to reduce 
all integrals (15. lip , with arbitrary integer powers of denominators and arbitrary (non- 
negative) powers of the selected scalar product in the numerators, to linear combinations 
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of 6 master integrals 




(5.12) 



This algorithm has been implemented in the SCHOONSCHIP [20] package Mincer [21] and 
later re- implemented [22] in FORM [23^1 . It has also been implemented in the REDUCE [2^ [25] 
package Slicer [26]. Only the last, non-planar, topology in (15. lip involves the last, non- 
planar, master integral in (I5.12p . 

The first 4 master integrals are trivial: they are expressed via Gn (15. 5p . and hence via 
F-functions. The 4-th one differs from the 3-rd one (G3) by replacing the two-loop subdi- 
agram: the second one in (15. 4p {G2/ {—k'^Y^'^) by the first one {G\/ {—k'^Y~'^). Therefore, 
it can be obtained from G^ by multiplying by 

GlG{lA-d) _2d-5Gl 
G2G{l,3~d) ~ d-3G^' 

The 5-th master integral is proportional to the two-loop diagram G(l, 1, 1, 1, n) with a non- 
integer index of the middle line n = 2 — d/2. The 6-th one, non-planar, is truly three- loop 
and most difficult. 

6 HQET self-energy diagrams 



There are 2 generic topologies of 2-loop HQET self-energy diagrams: 




(6.1) 



"'Unfortunately, Mincer does not produce linear combinations of 6 master integrals (|5.12|) : recursively 
1-loop integrals are expressed via F functions and expanded in e = 2 — d/2, so that the contributions from 
the first 4 master integrals cannot be separated. 
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The denominators of the second one are hnearly dependent, and there is one irreducible 
scalar product in the numerator. All these integrals, with any powers of denominators 
(and with any power of the numerator of the second diagram), can be reduced [15] (see 
also [27]) to 2 master integrals 



where 




is the n-loop HQET sunset integral. 

There are 10 generic topologies of 3-loop HQET self-energy diagrams: 




Diagrams in the first two rows have one scalar product which cannot be expressed via 
denominators; those in the third row have one linear relation among heavy denominators, 
and hence two independent scalar products in the numerator; those in the last row have 
two relations among heavy denominators, and hence three independent scalar products 
in the numerator. All these integrals, with any powers of denominators and irreducible 
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numerators, can be reduced [27] to 8 master integrals 




using integration by parts. This reduction algorithm has been implemented as a REDUCE 
package Grinder [27]. The first 5 master integrals can be easily expressed via F functions, 
exactly in d dimensions. The next 2 ones reduce to two-loop ones with a single (i-dependent 
index; the last one is truly three- loop. 

7 Equivalence of IBP relations for integrals with the 
same total number of loop and external momenta 

IBP relations for integrals with the same M = L + E but different E can be made equiva- 
lent [28] (then they differ only by boundary conditions — sets of trivial sectors). In other 
words, the IBP relations for an L-loop Feynman integral fl2.8l) with E external momenta 
can be reduced (after some re-definitions) to exactly the same form f l3.17p as for the M-loop 
vacuum integral (in which we also integrate in df^pi ■ ■ ■df'pE)- 

Consider some L-loop integral /(rii, . . . ^tln) with E external momenta at some kine- 
matic point Sij = > L, j > i). We number all such (i,j) pairs by integers 
a G [ + 1 , ] in some way, where 

K^N + N.^'J^ ,7^1) 

is the number of all scalar products Sij of M vectors Qi {j > i), and introduce the quantities 
Da = —Sij + s'^j = —Sij+ml (a G [A^+1, K]) for these pairs. Now all the quantities Da 
can be written in a uniform way (12. 6p (with both sums up to M). All the scalar products 
Sij can be expressed via Da in a uniform way (12.71) with the sum up to K. 

Let's expand this integral /(ni, . . . ,n]y) (12.81) in a formal series in Sij — {L < i < 
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J < M): 

oo oo 

/(ni, . . . ,nAr) = ^ ■ ■ ■ /(rii, . . . , wat, hn+i, • • • , nK)£'^^y~^ ■ ■ ■ D'^^~^ . (7.2) 

njv+i=l nK=i 

In fact, we are mainly interested in the value of our integral at this kinematic point {un+i = 
■ ■ ■ = uk = 1), the derivatives (given by higher values of these indices) are not our primary 
goal. We assume J = if any index tin+i, . . . , is < 0, so that there is only one sector 
with respect to each of these indices. 

Applying operators Oij with z G [1, L], j > i to the integrand /, we get the usual IBP 
relations (13.171) . In other words, the "evolution" of /(ni, . . . ,71^,11^+1, . . . ^uk) with re- 
spect to the "internal" indices rii, . . . , n^v is governed by "vacuum" IBP relations (13.171) . 

We need Ne additional relations which govern the "evolution" with respect to the 
"external" indices n^v+i, . . . , tik- To derive them, we apply A^; operators Oij = qj ■ di 
{L < i < j < M) to the definition (17. 2p . In the left-hand side, we get 

N N L 

^TAl.k-' (b- - ml) I{n,, . . . , n^) ; 

a=l fe=l m=l 

here each integral can be expanded like (17.21) . In the right-hand side, Oij acts on the 
product of Da with a E [N + 1, K]. To get the coefficient of -Djv+Y ^ ' ' ' w6 shift 

the index — )■ + 1 in terms with I /Da, and — )■ — 1 in terms with Df,, and obtain 

K K M 

J2 Yl XI ^^^^^i(^ + '^™Kb~ -"^6)a+/(ni,...,n^,nAr+i,...,n^) 

a=N+l b=N+l m=L+l 

(note the opposite order of the operators). In order to combine this with the left-hand 
side, we commute the operators (I3.15P and use 



Finally, the relations for the "external" indices are 

K K M 

E + 1)5,, + (b- - ml 

a=l b=l m=l 

X I{ni, . . . , riAT, riN+i, • • • , uk) = . 



(7.3) 



They are similar to the "internal" ones (I3.17p . but contain E + 1 instead of d. 
We can make these relation exactly "vacuum" ones using the substitution 

I{n,,...,nN) = D'^'^+'-'''>/'l{m,...,n^), 

^ det Sj,- . . ,^ , , (7.4) 

D = -nrJ^ z,jG[L + l,M] 



det4 
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The new functions I{ni, . . . ,nN) are expanded in Da with a > in the same way as 
in (17. 2p . with the coefficients I{ni, ... ,11^,71^+1, ... ,nK)- At Sij = s°- / coincides with 
/: i{ni, . . . , nAT, 1, . . . , 1) = /(rii, . . . , nA?, 1, . . . , 1); higher /(ni, ...,nN, un+i, • • • , rix) are 
hnear combinations of I{ni, . . . , n^r, 't-at+i, • • • , n^). The operator Ojj can be written as 

M Q 

m=L+l 

The derivatives of are 

d,jD = (2 - ■ (s-i),, G [L + 1, M]) 

(the extra factor 2 at z 7^ j comes from the fact that Sij appears in the determinant D 
twice, at the positions ij and ji). Therefore, 



= {E + 1- d)D^^+^-''^'H, 



I] ; 



and an additional term {d — E — l)Sij appears in (17.31) . All relations for 
J(ni, . . . , n^v, njy+i, . . . , uk) have the same form as for the M = L + E loop vacuum integral 
Io{ni,...,nK). 

The boundary conditions may differ. In the case of /(ni, . . . ,11^,11^+1, ■ ■ ■ ,nK), all 
sectors with Ua < (a G [N + 1,K]) are trivial; this does not have to be so for the 
corresponding vacuum diagram. 

Now we shall consider a few examples. 

Let's consider the 1-loop self-energy diagram (Fig. [TTk : we set m = 1): 

^(^i'^2) = -^y"-^0^, D, = l-{k + p)\ D2 = -e. (7.5) 

This integral vanishes at ni < 0. Suppose we want to calculate it on the mass shell = 1. 
We introduce D3 = 1 — p^, and re-express the integral as 

M(ni, n^) = (p^)(2-'^)/2^(ni, ns) , (7.6) 

according to (17.4p . Then we expand it in D^: 

00 

M(ni, ri2) = ^(^1' ^2, ^3)^3'"' (7.7) 

713 = 1 

(M(?7.i, ?T,2, ""-3) vanishes at < 0). There is 1 master integral: 

M(ni, n2, ris) = c(ni, na, n3)M(l, 0, 1) . (7.8) 

The d ■ k and d ■ {k + p) IBP relations for M(ni, 722) are 

[rf - m - 2n2 + nil+ {D-s - 2")] M(ni, ^2) = , 

[d - 2ni - n2 + 2nil+ + n22+ (D3 - 1")] M(ni, ^2) = 
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Figure 11: Massive 1-loop self-energy and 2-loop vacuum diagram 



(here is an external kinematic parameter). For the expansion coefficients this means 

(7.9) 



[c? - ni - 2n2 + nil+ (3" - 2")] M(ni, ^2, ng) = , 
[d - 2ni - n2 + 2nil+ + n22+ (S" - l")] M(ni, ^2, n^) = 



(we shift the summation index 77,3 —t- 723 — 1 in the terms with D^). In order to find the 
"evolution" in n-^, we apply p ■ d/dp to 07.71) : 

p ■ -^M{ni, n2) = [-ni + nil+ (2 + 2" - D3)] M(ni, ns) , 
op 



ni, 722,^3 



"3 = 1 



)p-^D-r' = 2 J2 M(ni,n2,n3)(n3 - 1) {0^-' - ^r^') 

"3 = 1 



2 ^ [(n3 - 1 - n33+) M(ni, n^, n^) 

"3 = 1 



r)"3-i 
^3 



Collecting these pieces together, we get 

[d-m- 2n3 + nil+ (2 + 2- - 3') + 2n33+] M(ni, n2, ^3) = . 
Now let's consider the 2-loop vacuum diagram (Fig. [TTb . m = 1): 

1 /• d%d'^k2 



(7.10) 



V"(ni,n2,n3) 
Di = 1 - (A;i + 



(7.11) 



^2 • 



This integral vanishes at rii < or < 0. The di ■ ki, di ■ {ki — ^2)? and 82 ■ k2 IBP 
relations for y(ni, ?t,2, 'n.3) have exactly the same form as 07.91) . 07. 101) . as expected. There 
is 1 master integral: 

V{ni,n2,n3) = co(ni, n2, n3)V(l, 0, 1) . (7.12) 
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In this problem, the boundary conditions for M{ni,n2,ns) and V{ni,n2,n3) coincide, 
and hence c(ni, n2, ris) = co(ni, n2, na). Therefore, the 1-loop on-shell self-energy diagram 
M{ni,n2) = M{ni,n2, 1) is related to the 2-loop vacuum diagram V{ni,n2, 1) (in which 
the index of the "former external" line is ns = 1): 



M{ni,n2) _ V{ni,n2, 1) 
M(1,0) " V^(1,0,1) 



(7.13) 



Explicit expressions for M(ni, 712) and V{ni,n2, n^) can be found, e. g., in the textbook [5]. 
It is easy to check the relation f l7.13p using 



T(z)T(l-z) 



TT 



smvr^; 



Reduction of 2-loop on-shell self-energy diagrams to master integrals has been consid- 
ered in I29]: 



— {^^^^^^1^^^ = Ci(ni, . . . ,725) 



+ 02(^1, . . . ,^5) 




(7.14) 




C3(ni, . . . , ns) -^ie- + C4(ni, ...,215) 




+ C5(ni, . . . ,^5) 




(7.15) 



These IBP reduction relations are equivalent to those for 3-loop vacuum diagrams (also 
considered in |22]): 




q(ni, . . . ,^0) 




+ ^(ni, . . . ,^6) 




+ C2(ni, ... ,726) 




+ 04(^1, . . . ,n6) 



(7.16) 




(7.17) 
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However, boundary conditions differ: on-shell integrals vanish at ne < 0, but vacuum 
integrals don't. Therefore, it is not possible to obtain Cfe(ni, . . . , ns) from c'^(^i, . . . , n^-, 1). 
Reduction for both classes of diagrams is implemented in a REDUCE package Recur sor |29j . 
For 2-loop on-shell diagrams, ng = 1 is specified, and a special control flag is set which is 
responsible for nullifying all integrals with ne < during reductiori^l. 




Figure 12: 1-loop vertex diagram (all lines are massless) 

Our next example is the 1-loop massless vertex diagram (Fig. [T2l) with 2 on-shell legs: 
Pi = P2 = ^- Its IBP relations are equivalent to those for the 2-loop massless self-energy 
diagram (Fig. |8]), because both are equivalent to the same 3-loop vacuum diagram with 
1 massive line. These 2-loop self-energy diagrams are expressed via 2 master integrals 
(Sect. ED : 




(7.18) 



The boundary conditions for the 1-loop vertex are different: integrals with n4 5 < vanish. 
Such integrals produce the second master integral in fIT.lSp . This leads to a very simple 
prescription: replace the second master integral in f l7.18p by and the flrst one by the 
vertex master integral. 




(7.19) 



^Thc redefinition (17.41) is not used in this package; therefore, some terms in the vacuum IBP relations 
are absent in the on-sheh ones. Their inclusion is controlled by the same flag. 
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Similarly, reduction for 2-loop massless vertex diagrams with pi = p 2 = is equivalent 
to 3-loop massless self-energy diagrams [28]. This method is applicable to diagrams with 
any number of external legs. 

8 Grobner-bases methods 

For complicated problems, constructing reduction algorithms by hand becomes impracti- 
cal, and it is desirable to automate this process. Some of such systematic approaches are 
based on Grobner bases. A simple introduction to Grobner bases for systems of polyno- 
mial equations with commuting variables is given in Appendix [Bl In the IBP reduction 
problem, polynomials of non-commuting operators are used. An approach using the shift 
operators f l3.10p was proposed in |30| |3T] and implemented in a Mathematica program 
FIRE [32] . 

Suppose we are considering the sector rii < 0, n2 > in Fig. [131 Any integral in the 
sector can be expressed via the integral at its corner using the shift operators 2+: 

/K,n2) = (l-)-"^(2+)"-V(0,l). (8.1) 





ni 



Figure 13: Normal form of IBP relations in a sector 

All IBP relations in this sector can be written in the normal form containing only 1~ 
and 2+: we act on the original form of these relations by 1~ and 2+ sufficiently many times 
to get rid of l"*", 2~ (but not too many times, see Fig. [13]). They are 

Q,,,(nJ(l-)^'^(2+y^~0 (8.2) 

jl J2>0 
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(this means that these operator polynomials give when applied to J; note that ria don't 
commute with the shift operators). We need to fix some total order of the monomials 
constructed from 1~ and 2"'"; this is equivalent to fixing a shift-invariant total order for 
the integrals in our sector (Sect. [3]). Now we can construct the Grobner basis for the 
system (18. 2p (more exactly S'-basis [301 Elj), reduce the monomial (1~)~"^ (2+)"^~^ (18. ip 
with respect to it, and apply it to /(0, 1). Integrals from lower sectors are considered trivial 
and already known. Irreducible monomials give master integrals in this sector. 

The idea to use Grobner bases for IBP reduction was originally proposed in [33] in 
a somewhat different setting. Let's assume that each line has a separate mass ma, and 
consider integrals without numerators (numerators can be eliminated by shifting d [5]). 
Differentiation in is equivalent to the raising operator (I3.13P : 

^ -a+ . (8.3) 



dm 



a 



Any integral in the positive sector can be expressed via the integral at its corner: 



/K,...,„0 = (-a;^j ...^-_) (8,4) 



The IBP relations can be written as 



J]C,,..,Jm?,...,m^)(^J ■■■[g^] -0. (8.5) 



Defining a total order for monomials constructed from d/dm\, we can find the Grobner 
basis of the system (18. 5p . Then we can reduce the monomial in (18. 4p with respect to this 
basis. Irreducible monomials give the master integrals. 

This method has been implemented in Maple and successfully applied [34] to the prob- 
lem of reduction of 2-loop self-energy diagrams with all 5 masses different and arbitrary 

[35] . If there are many zero (or equal) masses in the problem, this approach requires to 
solve a more difficult problem with all masses being different first, and this may lead to 
very lengthy intermediate expressions. 

9 Baikov's method 

The Feynman integral (12. 8p can be written as an integral in scalar products SjJ§. The 
integration measure is 

dVk2 ■ ■ ■ d'^kL = d'^-^k^\d''-'^+^k^^d^'-^k2\\d''-'''+^k2^ ■ ■ ■ d^'-'^kLlld'^-^'+'^kLl. , (9.1) 

where ki\\ lies in the subspace spanned by /c2, • • • , ^l, Pi, ■ ■ ■ , Pe', k2\\ lies in the subspace 
spanned by fcs, . . . , ki, pi, . . . , pe; and so on. 



^Here we consider the Euclidean case, in order not to have comphcations with signs and ±iO. 
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The volume elements are 



tM-it. _ dsi2dsi^ ■ ■ ■ dsiM 

d'''-''koU = 



,M-2r dS23dS24:- ■ ■ dS2M 



Gy^{h,...,kL,pi,...,PE)' (9.2) 

jM-Lt, _ d'^L,L+ldSL,L+2 ' ' ' dSiM 

where 

G{qi, ...,qn) = det qi ■ qj (9.3) 

is the Gram determinant {G^^^{qi, . . . , g„) is the volume of the parallelogram formed by qi, 
qn)- 

We can integrate over the angles in d'^ki±: 

d ki± = —^nk^j_ dk^j_ , 

where 

27r"/2 

is the n- dimensional full solid angle. Then we replace dk'fj_ = dsu; ki± is the height of the 
parallelogram with the base formed by fcj+i, . . . , ki, pi, . . . , pe and the extra vector ki, and 
this is the volume of the whole parallelogram divided by the area of its base. Therefore, 

,d-M+li, _ (G{ki,...,kL,Pl:...,PE)\~'^^'^^^'^,^ 

2 \G{k2,...,kL,Pi,...,PE)J 

,d-M+2i, _ ( ^(^2, • • • , fcL,Pl, • • • ,Pg) V"' ^^^''^ 

2 \G{ks,...,kL,Pi,...,PE)J (9.5) 

d'^-^-'k,, = W,,,, ( G{k,.Pu....PE) \ ^^^^ _ 



G{pi, ...,Pe) 

All the Gram determinants except the first and the last ones cancel in the measure (19. ip . 
and [6] 

I{ni, . . . , n^v) = j d'^ki ■ ■ ■ d'^kif = 

„~L(L-l)/4-L£;/2 

-G(p„...,P^)(-^+^+^)/^ 



nLr(-f±i) ^-'•••'-^ (9.6) 
X / l[l[ds,,G{ku...,kL,Pu...,PEY''"''-'^^'f. 

i=l j=i 
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The integration region has a comphcated shape; the Gram determinant vanishes on its 
boundaries. This formula can be used as a definition of the c?- dimensional integral. Usually, 
another definition based on the a (or Feynman) parametrization is used. But then many 
simple properties, like the possibility to cancel identical brackets in the numerator and the 
denominator of the integrand /, become theorems needing non-trivial proofs. Here the 
possibility to cancel brackets (depending on Sjj) is obvious. 

The integration variables Xa = Da (a G [1, A^]) [36j can be used instead of Sif 

...,n^) = C I "^lY'^^l^ Pix, - m?, . . . , - m^)(^-^^-i)/2 , (9.7) 



X-\ ' ' ' oc 



N 



where the Baikov polynomial is 



N 



a=l 



and the normalization constant is 



_-L(L-l)/4-L£;/2 

(in the last determinant, the pair (i, j) with j > i is considered as a single index). 

Acting by the operator a~ on (19. 7p multiplies the integrand by x^; acting by a"*" replaces 
by na/x"^"-^^ = —dail/Xa"^), we integrate by parts and take into account the fact that 
the polynomial P vanishes at the boundaries: 



a I = C 



f dXi ■ ■ ■ dXN p(d-A/-l)/2 
/ -r"! . . . -r"'^ 

-T_n f '^^1 • • • d^N ^ p(d-M-i)/2 

J Xi ■ ■ ■ Xjq 



(9.9) 



Note that [5a, x^] = 5abi in accordance with fl3.15p . The IBP relations Ojj(a+,a )/ = 
become ES] 

0,,{da.Xa)P^'-''-'^"{xa - ml) = Q . (9.10) 
Let's check that P{xa) (19. 8p indeed satisfies this requirement. If j < L (g^ = kj), 



N N M 



0,,(9a, Xa) = d6,, ^r^^,(l + S^^)^a{xt - ml) 

a=l b=l m=l 

N N M 

= ^d-M- 1)6., - E E E ^r<(l + Srm){Xb - ml)da 



a=l 6=1 m=l 

(see (I3.17p . note the order of operators). We substitute 

N Q N 

"y^, v4™9a = dmi = ; ^mji^b — m^) = Smj 

a=l 6=1 



30 



and get 

M 

Oij = {d- M - l)5ij - + 6mi)Smjdmi 

m=l 

(this form can be directly obtained from fl3.7p if one takes into account that the order of 
Smj and dmi in fl3.7p is interchanged due to integration by parts with respect to Xa in fl9.9p ). 
The derivatives of P = det s,,- are proportional to elements of the inverse matrix s~^: 



d,jP = (2 - 6,j)P ■ {s 



And now we see that fl9.10p is satisfied. The proof at j > L is similar, the only difference 
is that Sij with i,j>L (external kinematic quantities) are not expressed via Xa but kept 
intact. 

The Feynman integral (12. 8p can be expressed as 

I{n) = ^Ck{n)I{nk) , Ck{ni) = 6m , (9.11) 

k 

where are the indices of the master integrals. The coefficients Ck{n) obey the same IBP 
relations 

0,,(a+,a-)cfc(n) = (9.12) 

as the integral I{n), but with different boundary conditions. Let's consider a sector 
'^i,...,m > 0, nm+i,...,N < (any sector has this form after a suitable re-numbering of ra's). 
Suppose there is 1 master integral in this sector: its corner = 1, nm+i,...,N = 0. Then 

only integrals from this sector and higher ones can contain this master integral I{nk). In 
other words, the coefficient Ck{n) vanishes if any index ni ... ^ is < 0. 



Let's consider [271 

dx\ f dxui f dxyyi-\-i^ f dx j\j 



p(d-M-l)/2^ (9.13) 



where the contours in the first m integrals are small circles around the origin. This integral 
satisfies the IBP relations due to fl9.10p (because boundary terms from integration by parts 
in Xa vanish), and vanishes when any index ni ... ^ is < 0. It is natural to assume that it 
is a linear combination of Cfc(n) and Ck'{n) for master integrals in lower sectors (they also 
vanish when any of the indices rii ... ,„ is < 0). In other words, Cfc(n) is a linear combination 
of fl9.13p and similar integrals where some more / are replaced by ^. Coefficients in this 
linear combination are fixed by the boundary conditions (19. lip . 

As a simplest example, let's consider the 1-loop massive vacuum diagram (Fig. H]) (with 
m = 1). In Euclidean space, D = k"^ + 1. The Baikov polynomial is P{x) = x, and 



V{n) 



C r^(^-l)^'^-'^/' (9.14) 
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(it is easy to check that this formula reproduces the standard result for V{n)). There is 1 
master integral V(l): V{n) = c(?2)V"(l). The coefficient c(n) is 



c n 



dx 
1 



[X 



{n-l)\ 



\{d-2)/2 



d 



(n — 1)! \dx 



d 

V n 

2 



ri-l 



[X 



\{d-2)/2 



x=0 



(9.15) 



Taking into account c(l) = 1, we reproduce fl3.2ip . 

Now let's consider the 1-loop massless self-energy (Fig. [7]). In Euchdean space, Di = k^, 
D2 = {k + (we set = 1), therefore the Baikov polynomial is 



P{Xi,X2) 



Xi 

X2-XI—1 
2 



2 
1 



(X2 -Xi-1] 
4 



(9.16) 



There is only 1 non-trivial sector, and 1 master integral G{1,1). The coefficients of this 
master integral are 



c(ni,n2 



dxi f dx 



X 



n2 



P(Xi,X2 



Ad-3)/2 



(ni - 1)! (n2 - 1)! \dxi 



d 



ni— 1 



dX'} 



P(Xi,X2 



Ad-3)/2 



(9.17) 



xi=X2=0 



Application of this formalism to reduction of 3-loop vacuum diagrams is discussed in 
detail in [3S]- It was used in the formidable task of reduction of 4- loop massless self-energy 
diagrams [TU]. Some additional examples can be found in [32] and Chapter 6 of [2]. 

Using this formalism, P. A. Baikov [IQ] has constructed an elegant proof that the 3- 
loop massless non-planar integral (the last one in (16. 5p ) cannot be reduced to lower sectors. 
Indeed, suppose such a reduction exists, i. e. the equation 




Cl 





(9.18) 



is a linear combination of the IBP identities. Then it must hold not only for the Feynman 
integrals I{ni, . . . ,ns,ng), but also for any solution of the IBP relations (the indices ni § 
correspond to 8 denominators, and ng to the irreducible numerator). Let's apply it to 



s(ni, 



f dxi f dxs f dxg 



(d-5)/2 



This integral vanishes if any of the indices ni . g is < (all terms in fl9.18p except the first 
one are absent). For the first term we have 



s(l,...,l,n9) 



^-^Pio,...,o,xgr-'y^. 
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where P(0, . . . , 0, xg) ~ Xg(l — xg)^ (for a suitable choice of the irreducible numerator 
Dg [IQ]). This means that we may choose the integration contour in xg to go from 
to 1 (so that there are no boundary terms, and s satisfies the IBP relations). Then 
s(l, ...,1,0) 7^ 0, and we have a contradiction. This means that flQ.lSp cannot follow from 
the IBP relations. 

Why doesn't a similar reasoning apply to the planar diagram (the second one in (16. 4p ) 
which is reducible? For this diagram, P(0, . . . , 0, Xg) ~ Xg, and we cannot choose a suitable 
integration contour for Xg. 

A form of the irreducibility criterion suitable for application in complicated problems 
(such as 4- loop massless self-energies) has been proposed in [H]. 

10 Conclusion 

Approaches to the problem of reduction of Feynman integrals can be classified as following: 

• Generic 

• Construct an algorithm and implement by hand: Mincer, . . . 

• More automated approaches 

• Grobner-bases methods 

• Lie-algebra based method 

• Baikov's method 

• Specific numeric Tii. Laporta algorithm (AIR, FIRE, Reduze. . . ) 

The most straightforward approach to reduction is to substitute specific integer values 
for all the indices Ua and to solve the resulting huge linear system. If we consider some 
region of size R around the origin, there is L{L + E) relations per point in this region, 
and only 1 unknown per point (plus some unknowns outside the region near its surface, 
but their number is proportional to the area of the surface, and becomes negligible as 
compared to the volume of the region at sufficiently large R). The system of IBP relations 
is highly redundant. In each sector, we can use these relations to reduce more complicated 
integrals to simpler ones (with respect to some total order); integrals from lower sectors 
are considered trivial and already known. A few integrals which cannot be reduced any 
further are the master integrals. 

This is called the Laporta algorithm [12]. It may require solving huge linear systems. 
There are several publicly available implementations: AIR |33] (in Maple), FIRE [32] (in 
Mathematica), Reduze [H] (in C++, using the GiNaC library [35]). There are also nu- 
merous implementations which are not publicly available (several in Karlsruhe, several in 
Edmonton, by M. Czakon, etc.; FIRE version 2 is rewritten in C++ (with Fermat [l6]), but 
it is also not public). 

Many implementations of the Laporta algorithm are written in C++; for algebraic op- 
erations they either use GiNaC, or talk to an external CAS (e.g., Fermat). They have to 
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use very large number of expressions; in order not to exceed the size of the main memory, 
key-value databases (e.g., Kyoto Cabinet [37j) may be used. 

It is possible to reduce the redundancy of the system of IBP relations [14] using its Lie 
algebra structure (13. 6p . The set of relations 

di ■ ki+i {i e kL+i = ki) ; 

(10.1) 

y^^dj-kj 

i=l 

is sufficient: all the other Oij can be obtained by commutators from this smaller subset. It 
is still redundant: there are L + E + 1 relations per point, but this is smaller than L{L + E) 
relations per point in the complete set. 

The Laporta algorithm is universal. However, it requires one to solve very large linear 
systems, if the indices \na\ are not small, because the dimensionality of our integer space 
is large. Also, if only one \na\ ~ i? is rather large, the algorithm typically requires one 
to consider a region of size ~ i? in all directions around the origin. If a special-purpose 
algorithm for a family of Feynman integrals can be constructed, then reduction can be 
done much more efficiently, and higher \na\ are attainable. Also, reduction of integrals 
with one large \na\ ~ i? does not require to consider a huge number of integrals with all 
\na\ ~ R. 

Historically, reduction algorithms were constructed by hand, and implemented in vari- 
ous computer algebra systems (CASs). The pioneering program of this kind is Mincer [211 
[22] . It was used (and is still being used) for solving many physical problems. Some exam- 
ples of reduction programs for various classes of Feynman integrals are listed in Table [H 
It is impossible to produce a complete list, because special-purpose reduction programs 
for various diagrams were written and used in very many papers; the Table just contains 
typical examples. 

In recent years, the problems being considered became much more complicated (4-loop 
vacuum and self-energy diagrams). They cannot be solved using this traditional approach. 
Several attempts to automate construction of IBP reduction algorithms were made (see 
Sects. El El E])- Unfortunately, the problem has not been completely solved. It seems to 
be a well-defined mathematical problem, and I do hope that some universal and elegant 
solution will appear. 

A short comment on CASs used for implementing IBP reduction algorithms is in order. 
They form a wide spectrum. On one side, there is Mathematica: very convenient, with 
huge amount of built-in mathematical knowledge, with an advanced GUI; but very hungry 
with respect to both memory and CPU time (and also expensive). On the opposite end 
of this spectrum there is FORM [23j: low-level, nearly no built-in mathematical knowledge, 
but efficient and suitable for huge calculations. It has been recently released under GPL. 
REDUCE [2^ |25] is in the middle: a convenient high-level language, a lot of mathemati- 
cal knowledge (integrals, expansion in series, and much more), much more efficient than 
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Table 1: Some programs for reducing specific classes of Feynman integrals 



Program 


Ref. 


CAS 


Diagrams 


Mincer 


ED 

[22] 


SCHOONSCHIP 
FORM 


3-loop massless self-energies 


Slicer 


[26] 


REDUCE 


Recursor 


[29] 


REDUCE 


2- loop massive on-shell self-energies 

3- loop massive vacuum diagrams 


SHELL2 
0NSHELL2 


m 
m 


FORM 


2-loop massive on-shell self-energies 


MATAD 


m 

[52] 


FORM 


3-loop massive vacuum diagrams 


Grinder 


[27J 


REDUCE 


3-loop HQET self-energies 


SHELLS 


m 


FORM 


3-loop massive on-shell self-energies 


Tarcer 


[35] 

m 


FORM 

Mathematica 


2-loop self-energies with arbitrary masses 


Loops 


m 


REDUCE 


2-loop massless self-energies 



Mathematica ; it has no fancy GUI. A few years ago it became free software (BSD license). 
If the problem being considered is so large that the size of expressions is larger than the 
main memory of the computer, FORM is the only choice. It can work with expressions stored 
on disk efficiently. If any other system begins to swap, the situation is hopeless. On the 
other hand, if the problem fits in the main memory, REDUCE is a very reasonable candidate. 
Its language is well suited for this kind of problems (I've written a large package Grinder 
in it). 

Several years ago, an interesting benchmark was run on many different CASs [SB]. It 
was multiplication of large sparse multivariate polynomials. Unsurprisingly, specialized 
polynomial systems (Pari ^48j, Format [l6]) were at the top. The fastest general-purpose 
system was REDUCE (closely followed by maxima compiled with CMUCL, an efficient com- 
mon lisp compiler). FORM was 4 times slower than REDUCE; Maple 10 times slower; and 
Mathematica nearly 20 times slower. REDUCE has good implementations of polynomial 
GCD (important for working with rational functions, and hence for IBP reduction) and 
of polynomial factorization, including multivariate (both of these things are impossible or 
very difficult in FORM). REDUCE was the most widely used CAS in physics at the dawn of 
the CAS era, but its use has dramatically reduced later; now, when it's free software, it 
would be good to use this excellent system more often. But, as I said, FORM has a unique 
advantage: it is the only system that can work with expressions which are larger than the 
main memory of the computer. 

Acknowledgements. I am grateful to P. A. Baikov, K. G. Chetyrkin, A. V. Kotikov, 
R. N. Lee, A. V. Smirnov, and V. A. Smirnov for reading the manuscript and suggesting 
improvement; to D. J. Broadhurst and M. Steinhauser for discussions of various aspects of 
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A Dimensional regularization 

During these lectures, we are going to live in a c?- dimensional space-time (1 time and d — 1 
space dimensions, the signature H — ■■■—). The dimensionality ci is a symbol; comparing 
it to integers (rf > 4 or < 4) makes no sense. If you want to consider a 100- dimensional 
subspace of the d- dimensional space, this is perfectly legal; there will also be a (c? — 100)- 
dimensional subspace orthogonal to your one. 

An L-loop Feynman integral / (12. 8p has dimensionality of mass to the power Ld — n, 
where n is an integer dimensionality of the denominators. If / contains no dimensional 
parameter, we can write no result for such an integral except 0. For example. 




(A.l) 



and so on. 

If / contains 1 dimensional parameter, its power is given by dimensions counting: 




(A.2) 

lid) , 



where I{d) are dimensionless functions of 1 variable d. 

If / contains 2 dimensional parameters, it is a non-trivial function of their ratio (in 
addition to d): 




Integrals with larger numbers of scales non-trivially depend on all dimensionless ratios. 
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B Grobner bases 



Here we briefly introduce Grobner bases for systems of polynomial equations with com- 
muting variables. For more detail, see, e.g., the textbook [57] . 

B.l Statement of the problem 

Suppose we have n variables xi, . . . , a;„. They are not independent, but satisfy some 
polynomial equations pi = 0, . . . , = {pj are polynomials of Xi). Let's consider some 
polynomial q of the same variables. It is natural to ask if this polynomial is equal to due 
to the constraints on our variables or not. If there is another polynomial there is the 
question of their equality. 

These questions would become very easy if we had an algorithm reducing polynomials 
of dependent variables to a canonical form. Two equal polynomials reduce to the same 
canonical form; a polynomial equal to reduces to the canonical form 0. 

We can try to use the equations pj = for simplifying the polynomial q, i. e. for replacing 
its more complicated terms by combinations of simpler ones. But to do so we first have to 
accept some convention which terms are more complicated and which are more simple. 

B.2 Monomial orders 

We need a total order of monomials (i.e. products of powers of the variables x"^ ■ ■ -x"^"). 
An order is total if for any monomials s and t either s < t, or s > t, or s = t is true. An 
order is admissible if two properties are satisfied: 

• 1 < s for any monomial s; 

• if s < t then su < tu for any monomial u. 

One of the most popular monomial orders is lexicographic. Anybody who has ever 
seen a dictionary knows it. We are comparing two monomials: s and 

L — Xi X2 -^n ■ 

• nx> rrix^ s > t 

• rii < nil ^ s < t 

• Til = nil ^ 

• n2> ni2 =^ s > t 

• n2 < ni2 =^ s < t 

• n2 = ni2 ^ 
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Another popular order is by total degree than lexicographic. First we compare the total 
degree n = ni + n2 + - ■ ■+nn of the monomial s and the total degree m = mi + m2 + - ■ ■+mn 
of the monomial t: 

• n > m ^ s > t 

• n < m ^ s < t 

• n = m ^ compare lexicographically 
B.3 Reduction of polynomials 

Let's fix some admissible monomial order. We'll write polynomials in descending order: 
the leading term first, followed by the rest ones. We'll normalize all polynomials pj in such 
a way that the coefficient of the leading term is 1. Now they can be used as substitutions 
which replace the leading term by minus sum of the remaining ones. I. e. if some term of a 
polynomial q is divisible by the leading term of some polynomial Pi, we remove this leading 
term and insert minus sum of the remainder terms of Pi instead. This is called reduction 
of the polynomial q with respect to the set of polynomials pf, if none of the substitutions 
is applicable, the polynomial q is called reduced. 

For example, let's fix the lexicographic order with x > y, and consider the set of 
polynomials 

Pi = x^ + t/^-l, p2 = xy-\. (B.l) 

We want to reduce q = x'^y with respect to this set. This can be done in 2 ways: we can 
first reduce q with respect to pi or p2 (Fig. [H]). 




Figure 14: Reducing q with respect to pi and p2 

So, we have obtained two different results, qi and q2; both cannot be reduced further. 
In fact they are equal due to pi = and p2 = 0, but this is not evident. Every time when 
more than one substitution can be applied to a term of a polynomial q (in this particular 
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case, we can replace either or xy in x^y), a fork appears; maybe, its branches join later, 
but maybe, they don't (as in this case). This example motivates the following definition: 

A set of polynomials pi, . . . , p„ is called a Grobner basis (for a given monomial order) 
if reduction of any polynomial q with respect to this set is unique. 

This definition is not constructive: it does not say how to check if a given set of 
polynomials forms a Grobner basis. There exist equivalent constructive definitions. 

B.4 5- polynomials 

In our example, the constraints pi = and P2 = allow us to simplify the monomials 
x^ and xy. Do these constraints contain an extra information usable for simplification 
but not obvious? Yes, they do! Let's multiply pi and p2 by monomials (i.e. products of 
powers of variables) in such a way that their leading terms become identical (equal to the 
least common multiple of the leading terms of pi and p2)- Then we subtract the second 
polynomial from the first one. The leading terms cancel, and we get a new polynomial 
with a new leading term which can be used for simplifying terms in q (because this new 
polynomial also vanishes). This polynomial is called the S-polynomial >S'(j9i,p2) (from the 
word subtraction). In our example 

Pi = + — 1 = X ?/ x'^y + y^ — y = 
P2 = xy — ^ = XX x^y — |x = 

\x + y^ -y = Q 

This polynomial can be added to the system of constraints pi = 0, j92 = 0. Let's normalize 
its leading coefficient to 1: 

p^ = x + Ay^-Ay. (B.2) 

Now there is a new possibility for reduction (Fig. [T5l) . Reducing q2 with respect to ps, 
we again obtain qi. Reduction of q with respect to the set pi, p2, Ps leads to the unique 
result. In fact, pi, p2, Ps form a Grobner basis (though we have not proved this). 

The set of constraints pi = 0, p2 = 0, p^ = can be simplified by reducing these 
polynomials with respect to each other (and omitting identical zeros). Such a reduced 
Grobner basis is 

Pi = y^-y' + 4 = 0, P2 = x + Ay'-Ay = 0. (B.3) 

It has triangular structure: the first restriction involves only the junior variable y; the 
second one — y and x. This is a general feature in the case of lexicographical orders. 

B.5 Buchberger algorithm 

Generalizing this example, we can formulate an algorithm to construct a Grobner basis 
from a set of polynomials: 
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91 = 



-y^ + y 




Figure 15: Reducing q with respect to pi, p2, Ps 

• Given P = {pj} 

• S = {the set of pairs {pi,Pj) of these polynomials with i < j < n} 

• while S is not empty 

• choose and remove some pair {pi,Pj) from S; 

• calculate S'-polynomial S{pi,pj)] 

• reduce it with respect to P; 

• if the result is not 0, add this polynomial to P, 
and the corresponding pairs to S 

The set of pairs S alternatingly shrinks and grows. But it can be proved that this process 
terminates after a finite number of steps, and produces a Grobner basis P. 

C Using FIRE 

The Mathematica package FIRE [32] implements the Laporta algorithm and the Grobner- 
bases method (Sect. M)- It can be downloaded from the site [32]. You will need 2 files: 
FIRE. tar. gz and SBases_3 . 1 . .m (this version is newer than the one in the tar.gz file). 
Unpack the archive in some directory, copy SBases_3 . 1 . .m to it, and include this directory 
into your Mathematica $Path. 

Let's investigate the 2-loop massless self-energy diagram (Fig. |8]). The beginning of any 
program using FIRE is standard: 

«SBases_3. l.O.m 

«FIRE_3.4.0.in 

«IBP.m 



40 



Then we store the hsts of the loop and external momenta, and the list of the propagator 
denominators Da-, in global variables: 

Internal=-[kl,k2}; 
External=-[p}- ; 

Propagators={- (kl+p) ~2 , - (k2+p) ~2 , -kl~2 , -k2~2 , - (kl-k2) ~2} ; 
PreparelBP [] 

The set of Da must be linearly independent, and all the scalar products must be expressible 
via Da'i linear denominators are allowed. There is no special data type for vectors; kl"2 
and kl*k2 are interpreted as k\ and ki ■ k2- 

Let's have a look at the di ■ k2 IBP relation generated by the program: 

IBP [kl,k2] 

- a[3] + a[5] - p'^a[l]Y[l] - a[l]r [l]Fm[2] - a[l]r [l]Fm[3] - a[5]Y[5]Ym[3] 

- a[3]Y[3]Ym[A\ + a[5]Y[5]Ym[A] + a[l]r [l]rm[5] + a[3]r [3]rm[5] 

Here a[3] means n^, Y[l] means 1"^, and ym[2] means 2~. Individual IBP relations may 
contain scalar products of external momenta only linearly (here p^; if there were 2 external 
momenta pi and p2, the product piP2 could also appear). At this stage, their meaning is 
unambiguous. However, later, when combining many IBP relations, you can get, say, plp2, 
and you cannot know if it is plp2 or {pi ■ P2Y- Therefore, if is important to substitute 
some kinematic invariants for scalar products like piP2 immediately after generating the 
IBP relations. Here we list all the IBP relations: 

startinglist=-[IBP[kl,kl] ,IBP[kl,k2] ,IBP[kl,p] , 
IBP[k2,kl] ,IBP[k2,k2] , IBP [k2,p] }/ .p~2->-l ; 

(this list can be shortened to reduce the redundancy, see ( 110. ip ). 

The next step is to write down all the symmetries of the diagram. If the substitution 
of the line numbers 1, 2, . . . , — )■ si, S2, . . . , s^r is a symmetry, we include {si, S2, . . . , s^r} 
into the list. It is not sufficient to provide generators of the symmetry group, because the 
program does not try to multiply the elements of the symmetry group you provide; all 
elements, except the identity substitution, have to be included: 

SYMMETRIES=-[{2 ,1,4,3,5}, 
{3,4,1,2,5}, 
{4,3,2,1,5}}; 

Now we are ready to burn the fire: 

Prepare [AutoDetectRestrictions->True] 
Burn [] 
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These lines initialize data structures for using the Laporta algorithm. The option 
AutoDetectRestrictions->True instructs FIRE to detect trivial sectors automatically, 
using the Lee's criterion (Sect. [3]). It first appeared in SBases version 3.1.0; therefore, it is 
important to download this version. Alternatively, instead of autodetecting trivial sectors, 
you can provide their list by hand: 

RESTRICTI0NS={{-1 ,-1,0,0,0}, 
{-1,0,-1,0,0}, 
{-1,0,0,0,-1}, 
{0,-1,0,-1,0}, 
{0,-1,0,0,-1}, 
{0,0,-1,-1,0}, 
{0,0,-1,0,-1}, 
{0,0,0,-1,-1}}; 

Prepare [] 

Each element of the external list designates a set of sectors; —1 in the a-th position means 
sectors with Ua < 0; +1 means sectors with > 0; means "don't care", all sectors along 
Ua will do. Of course, in lists of trivial sectors (RESTRICTIONS) you will never use +1. It 
is important to include all trivial sectors; but if, say, ni < 0, ?7,2 < is sufficient for / = 0, 
you don't have to list 2^ sectors with respect to n3^...^5, just write in these positions. 

The integral /(ni, . . . , nTv) is denoted F [{nl , . . . ,nN}] in the program (don't forget that 
indices corresponding to irreducible numerators are always < 0). The result is computed 
in terms of the master integrals G[{nl, . . . ,nN}]. Let's calculate 

F[{1,1,1,1,1}] 

The program writes some messages explaining what's going on while applying the Laporta 
algorithm, and produces the result 

2(-10 + 3d){-8 + 3d)G[0, 1, 1, 0, 1] 2(-3 + d)G[l, 1, 1, 1, 0] 
(-4 + d)2 -4 + d ■ 

We see that there are 2 master integrals (15.41) . If things go well, you don't need to read 
messages detailing the progress of the algorithm; but they may be useful if something goes 
wrong. Let's try something more difficult: 

F[{2,2,2,2,2}] 

- (l8(-9 + d){-7 + d){-2 + d){-20 + 3d)(-16 + 3d){-U + 3rf)(-10 + 3rf)(-8 + 3d) 
G[0, 1, 1, 0, l])/((-10 + d)\-8 + rf)2(-4 + d)) 
4(-6 + rf)(-5 + d){-3 + d)(-708 + 242d - 27^2 + d^)G[l, 1, 1, 1, 0] 
^ {-10 + d){-8 + d) ■ 

FIRE can build S'-bases in some sectors, and later use them for reduction. If I add the 
lines 
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BuildBasis[-[l, 1, 1, 1, 1}] 
BuildBasis[{l, 1, 1, 1, -1}] 
BuildBasis[-[-l, 1, 1, 1, 1}] 
BuildBasis[-[-l, 1, 1, -1, 1}] 

between Prepare and Burn, then calculation of F [-[2,2,2,2,2}] takes 2 seconds instead 
of 20 (on my computer). 



Let's consider one more example: 3-loop vacuum diagram [29] in Fig. [T6l 



«SBases_3.1.0.ni 

«FIRE_3.4.0.ni 

«IBP.m 

Internal=-[kl , k2 , k3} ; 
External=-[} ; 

Propagators={l-kl"2,l-k2"2,l-k3~2,-(k2-k3)~2,-(k3-kl)~2,-(kl-k2)~2}; 
PreparelBP [] 

startinglist=-[IBP[kl,kl] ,IBP[kl,k2] ,IBP[kl,k3] , 

IBP [k2 , kl] , IBP [k2 , k2] , IBP [k2 , k3] , 

IBP [k3 , kl] , IBP [k3 , k2] , IBP [k3 , k3] } ; 
SYMMETRIES={ {2 ,3,1,5,6,4}, 




Figure 16: A 3-loop vacuum diagram 



{3,1,2,6,4,5}, 
{1,3,2,4,6,5}, 
{3,2,1,6,5,4}, 
{2,1,3,5,4,6}}; 



Prepare [AutoDetectRestrictions->True] 
Burn [] 

F[{1, 1,1, 1,1,1}] 



3(-10 + 3ci)(-8 + 3rf)G[l, 0, 1, 1, 0, 1] (-2 + dfG[l, 1, 1, 0, 0, 0] 
16(-4 + rf)(-7 + 2d) 8(-4 + rf)(-3 + d) 
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F[{2,2,2,2,2,2}] 



(27(-9 + d){-6 + rf)2(-4 + d){-22 + 3rf)(-20 + 3c;)(-16 + 3rf)(-14 + 3d) 

(-10 + 3d)(-8 + 3rf)(29 - 12rf + o, 1, 1, 0, 1]) 

/(l6384(-8 + d){-7 + d){-l7 + 2rf)(-15 + 2rf)(-13 + 2c/)(-ll + 2d)(-9 + 2rf)(-7 + 2d)) 
_ (-6 + d)\-A + d)^{-2 + rf)^(62 - Ihd + 1, 1, 0, 0, 0] 

512(-8 + d){-7 + c/)2(-5 + c/)(-3 + d) ' 

There are 2 master integrals fl7.16p . 
If I add the hnes 

BuildBasis [-[1,1,1,1,1,1}] 
BuildBasis [{1,1,1,-1,1,1}] 
BuildBasis [{-1 ,1,1,1,1,1}] 
BuildBasis [{1,1,1,-1,-1, 1}] 
BuildBasis [{-1 ,1,1,-1,1,1}] 
BuildBasis [{1,1,1,-1,-1 , -1}] 

between Prepare and Burn, calculation of F [{2 ,2,2,2,2, 2}] takes 30 seconds instead of 
100. 
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